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Abstract 

Given a finite planar graph, a grove is a spanning forest in which every component 
tree contains one or more of a specified set of vertices (called nodes) on the outer 
face. For the uniform measure on groves, we compute the probabilities of the different 
possible node connections in a grove. These probabilities only depend on boundary 
measurements of the graph and not on the actual graph structure, i.e., the probabilities 
can be expressed as functions of the pairwise electrical resistances between the nodes, 
or equivalently, as functions of the Dirichlet-to-Neumann operator (or response matrix) 
on the nodes. These formulae can be likened to generalizations (for spanning forests) 
of Cardy's percolation crossing probabilities, and generalize Kirchhoff's formula for 
the electrical resistance. Remarkably, when appropriately normalized, the connection 
probabilities are in fact integer-coefficient polynomials in the matrix entries, where the 
coefficients have a natural algebraic interpretation and can be computed combinato- 
rially. A similar phenomenon holds in the so-called double-dimer model: connection 
probabilities of boundary nodes are polynomial functions of certain boundary measure- 
ments, and as formal polynomials, they are specializations of the grove polynomials. 
Upon taking scaling limits, we show that the double-dimer connection probabilities 
coincide with those of the contour lines in the Gaussian free field with certain natural 
boundary conditions. These results have direct application to connection probabilities 
for multiple-strand SLE2, SLEg, and SLE4. 

1 Introduction 

1.1 Grove partitions 

A circular planar graph Q is a finite weighted planar graph with a set of vertices N on its 
outer face numbered 1, . . . , n in counterclockwise order. The vertices in N are called nodes, 
and the remaining vertices are called inner vertices. Define a grove to be a spanning 
acyclic subgraph (a forest) of Q such that each component tree contains at least one node. 
The weight of a grove is the product of the weights of the edges it contains. We study 
random groves where the probability of a grove is proportional to its weight. 

The term grove comes from Carroll and Speyer [CS04J and Petersen and Speyer [PS05J 
who studied a special case of (our) groves, with a particular family of underlying graphs. 
Since we are dealing with a natural generalization we chose to keep their terminology, and 
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refer to the special case they discuss as Carroll- Speyer groves, which we will discuss further 
in a subsequent paper [KW08] . 

The connected components of a grove partition the nodes into a planar (i.e. noncrossing) 
partition. For example, when n = 4, there are 14 planar partitions: 1234, 1|234, 2| 134, 
3)124, 4| 123, 12|34, 23|14, 1|2|34, 1|3|24, 1|4|23, 2|3|14, 2|4|13, 3|4|12, 1|2|3|4. There are no 
groves with the partition 13 1 24 because it is not planar (there is no way to connect nodes 1 
and 3 and nodes 2 and 4 by disjoint paths within a circular planar graph). For general n, 
the number of noncrossing partitions is the n th Catalan number C n = (2n)!/(n!(n + 1)!) (see 
[Sta99l ex. 6.19(pp), pg. 226]). 




7 . : 5 




Figure 1: A random grove (left) of a rectangular grid with 8 nodes on the outer face. In 
this grove there are 4 trees (each colored differently), and the partition of the nodes is 
{{1}, {2, 7, 8}, {3, 4, 5}, {6}}, which we write as 1 1 278 1 345 1 6, and illustrate schematically as 
shown on the right. 

If cr is a planar partition of 1, . . . , n, we let Pr(cr) denote the probability that a random 
grove of Q partitions the nodes according to o. Since groves with one component are trees, 
we refer to the partition a = 123 . . . n as the tree partition. When there are n nodes we 
call the partition a = l|2|3|---|n the uncrossing, since groves with this partition type 
contain no crossings (i.e. paths) connecting the nodes. We show how to compute for each 
planar partition a the probability Pr(cr) of it occurring in a random grove, as a function 
of the electrical properties of the graph Q, when Q is viewed as a resistor network with 
conductances equal to the edge weights. 

We let R it j denote the effective electrical resistance between nodes i and j, i.e., the 
voltage at node i which, when node j is held at volts, causes a unit current to flow through 
the circuit from node i to node j. (A good reference for basic electric circuit theory is 
[DS84].) Let Ljj denote the current that would flow into node j if node i were set to one 
volt and the remaining nodes set to zero volts. Though it is not obvious from this definition, 
Appendix A). The Lj,- are the negatives of the entries of the "response 



J .i-t 
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matrix" (Dirichlet-to-Neumann matrix) of (Q, N), on which we provide further background 
in |Appendix A[ see also |CdV98] . We prove 
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Theorem 1.1. For any planar partition a, 

Pr (a) integer- coefficient homogeneous polynomial in the Rij/2 's 
Pr(tree) where the degree is — 1 + Mparts of a, 

Pr(cr) integer- coefficient homogeneous polynomial in the Lij 's 



and 



Pr (uncrossing) where the degree is n — Mparts of a. 

This theorem is a corollary to our main theorem on groves, Theorem 1.21 in ISection 1.2[ 
which gives explicit formulas. These formulas are in effect a multinode generalization (for 
planar graphs) of Kirchhoff's formula Pr(l|2)/ Pr(12) = R\ 2 = \/L\ t 2 (which holds whether 
or not Q is planar) |Kir90j (see |Pem95l Thm 3.7] for a good exposition). (When n = 2 we 
have Li 2 = 1/Ri >2 , but this does not hold for n > 2.) We illustrate ITheorem 1.11 by giving 
the polynomials for n = 2 and n — 3, and in Appendix B we give the polynomials for n = 4 
nodes. For notational convenience we write 

Pr(er) := Pr(er)/ Pr(tree) and Pr (a) := Pr(cr)/ Pr (uncrossing). 

(The dots in Pr() remind us that in the normalization the nodes are disconnected.) The 
polynomials for n = 2 come from Kirchhoff's formula: 

Pr(12) = 1, Pr(l|2) = Ri, 2 , and Pr(12) = L 1>2 , Pr(l|2) = 1. 

When n = 3 the polynomials start to become more interesting: 

_ P7(123) = 1, P7(l|23) =±Ri, 3 + \Ri,2 - f#2,3, 

Pr(_2|13) = ii? 2)3 + \Ri,2 - |i2i, 3 , Pr(3|12) = \R^ + \R 2 ,z - \Ri, 2 , 

Pr(l|2|3) = ii?i i2 -Rl,3 + 2-^1,2-^2,3 + 2-^1,3-^2,3 ~ 4-^1,2 — 1^2,3 — 1^1,3' 

and 

Pr(123) = L12L13 + L lj2 -^2,3 + -^1,3-^2,35 
Pr(l|23) = L 2)3 , Pr(2|13) = Li >3 , Pr(3|12) = L lj2 , Pr(l|2|3) = 1. 

The n = 3 formulas also hold whether or not the graph is planar, but when n > 4 planarity 
becomes important. 

When n gets large these polynomials can have many terms, but for certain special classes 
of partitions, most notably the "parallel crossing" 1, n\2, n — 1(3, n — 2| • • • , the polynomials 
can be expressed in terms of determinants [CIM98] [FomOlj . We plan to discuss these and 
other determinant / Pfaffian formulas in a subsequent article [KW08J. 



1.2 The partition projection matrix and explicit formulas 

The formulas in ITheorem 1.11 in terms of the L^/s are the simplest ones to explain. For a 
partition r on 1, . . . , n we define 

F {i,j}&F 
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where the sum is over those spanning forests F of the complete graph on n vertices 1, . . . , n 
for which trees of F span the parts of r, and the product is over edges {i,j} of forest F. 

This definition makes sense whether or not the partition r is planar. For example, L11234 = 
-^2,3-^3,4 + £2,3-^2,4 + -^2,4-^3,4 and L 13 |24 = L 13 L 2i 4. As we shall see, the "L polynomials" of 
IThcorem 1.11 are in fact integer linear combinations of the L T 's: 

Pr((j) 



Pr(l|2| 



\n) 



(We write the superscript (t) to distinguish these coefficients from ones that arise in the 
double-dimer model in ISection 1.4[ ) The rows of the matrix V^' are indexed by planar 
partitions, and the columns are indexed by all partitions. In the case of n = 4 nodes, the 
matrix V® is 
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For example, the row for 1|234 tells us 
Pr(l|234) 



Pr(l|2|3|4) 
and the row for 12 1 34 tells us 

Pr(12|34) 



-^1|234 + -^13|24 — -^2,3-^3,4 + -^2,3-^2,4 + -^2,4-^3,4 + -^1,3-^2,4 



Pr(l|2|3|4) 



-^12|34 — -^13|24 — -^1,2-^3,4 — -^1,3-^2,4- 



We call this matrix "PW the projection matrix from partitions to planar partitions, 

since it can be interpreted as a map from the vector space whose basis vectors are indexed 
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by all partitions to the vector space whose basis vectors are indexed by planar partitions, 
and the map is the identity on planar partitions. For example, the column for 13 124 tells us 

13|24 projects to - 12|34 - 14|23 + 1|234 + 2|134 + 3|124 + 4| 123. 

(We could have written the right-hand-side as -ei 2 |34 - ei 4 |23 + ei|234 + e 2 |i34 + e 3 |i 24 + e 4 |i23, 
where the e's are basis vectors, but it is convenient to suppress the vector notation and 
instead write formal linear combinations of partitions.) 

The projection matrix may be computed using some simple combinatorial transforma- 
tions of partitions. Given a partition r, the r th column of may be computed by repeated 
application of the following transformation rule, until the resulting formal linear combination 
of partitions only involves planar partitions. The rule generalizes the transformation 






4—1 

13|24 -> -12|34 - 14|23 + 1|234 + 2|134 + 3|124 + 4|123 

(which is derived in ILemma 2.3I) to partitions r containing additional items and parts. If 
partition r is nonplanar, then there will exist items a < b < c < d such that a and c belong 
to one part, and b and d belong to another part. Arbitrarily subdivide the part containing 
a and c into two sets A and C such that a G A and c G C, and similarly subdivide the part 
containing b and d into B 3 b and D 3 d. Let the remaining parts of partition r (if any) be 
denoted by "rest." Then the transformation rule is 

AC\BD\ rest -> A\BCD\ rest + B\ACD\ rest + C\ABD\ rest + D\ABC\ rest 

- AB\CD\ rest - AD\BC\ rest . (Rule 1) 

In ISection 21 we prove 

Theorem 1.2. Any partition r may be transformed into a formal linear combination of 
planar partitions by repeated application of \Rule I\ and the resulting linear combination does 
not depend on the choices made when appluing \Rule 1\ so that we may write 



planar partitions a 



For any planar partition a, these same coefficients Va,r satisfy the equation 

Pv(a) 



Pr(l|2| 



\n 



partitions r 

for circular planar graphs. More generally, for any graph these coefficients satisfy 

Pv(t) 



t) 



Pr 1 2 ■■■ n) 

partitions t \ i i i / partitions r 

(The last equation specializes to the previous one because Vajr = 1, for planar r ^ a we 
have Va}- = 0, and for nonplanar r we have Pr(r) = for circular planar graphs.) 
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1.3 Double-dimer pairings 

An analogous computation can be done for the double-dimer model. Let Q be a finite edge- 
weighted bipartite planar graph, with a set N of 2n special vertices called nodes on the outer 
face, which we label 1, . . . , 2n in counterclockwise order along the outer face. A double-dimer 
configuration on {Q, N) is by definition a configuration of disjoint loops (i.e. simple cycles 
of length more than two), doubled edges, and simple paths on Q that connect all the nodes 



in pairs, which collectively cover the vertices of Q. See Figure 2 for an example. We weight 
each configuration by the product of its edge weights times 2 e , where i is the number of 
loops (a doubled edge does not count as a loop). The number of possible ways that the 2n 
nodes may be paired up with one another is the n th Catalan number (see [Sta99 > ex. 6.19(n), 
pg. 222]), and we show how to compute the probabilities of each of these pairings when a 
random double-dimer configuration is chosen according to these weights. We will show how 
to express these pairing probabilities in terms of polynomials in a set of variables that are 
analogous to the Lij variables for random groves. 




I 1 



4 



■ 



1 I 



3 1 



2 



i 



Figure 2: At the left is a double-dimer configuration on a rectangular grid with 8 nodes. In 
this configuration, the pairing of the nodes is {{1,8}, {3,4}, {5,2}, {7,6}}, which we write 
as 8 1 4 I 2 1 6 (odd-numbered nodes are always paired with even- numbered nodes). A double- 
dimer configuration is formed from the union of two dimer configurations, one on the graph 
qbw q g (dgfingd below and shown in the middle) and the other on the graph Q WB C Q 
(shown on the right). 



Let Z DD (Q, N) be the weighted sum of all double-dimer configurations. Since the graph Q 
is bipartite, we view the vertices as being colored black and white so that each edge of Q 
connects a black and white vertex, and we define the parity of a node to be the parity of 
its numerical label. Let Q BW be the subgraph of Q formed by deleting the nodes except the 
ones that are black and odd or white and even, and let Qff be defined as Q BW was, but with 



node i included in Qfj if and only if it was not included in Q BVJ , and similarly for node j. Let 
Z BW and ZfJ* be the weighted sum of dimer configurations of Q BW and Qf™ respectively, 



and define Z WB and Zf B similarly but with the roles of black and white reversed. Each 
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of these quantities can be computed via determinants, see |Kas67] and ISection 51 It turns 
out that Z DD = Z BW Z WB ; this is essentially Ciucu's graph factorization theorem |Ciu971 
Thm. 1.2] (except that in the factorization theorem, any edges that connect two nodes are 
reweighted by 1/2, and here they are not), see ISection 31 (The two dimer configurations in 
Figure 2] are on the graphs Q BW and Q WB .) The variables that play the role of Lj 7 - in groves 



are defined by 



X 



gBW j^B\N 



For each planar matching of the nodes a, let Pr (a) be the probability that a random 
double-dimer configuration has this set of connections. In ISection 3l we prove 

Theorem 1.3. For any planar pairing a on 2n nodes, 



Pr a) 



integer- coefficient homogeneous polynomial of degree n in the quantities Xij. 



This theorem is a corollary to our main theorem on double-dimer pairings, ITheorem 1.41 
in ISection 1.41 which gives the polynomials explicitly. To illustrate this theorem, we give a 
few examples. For notational simplicity let us define 



Pr(tr) = Pr(a)Z WB /Z 

In the simplest nontrivial case there are In = 4 nodes, and two possible planar pairings, 
{{1, 2}, {3, 4}} and {{1, 4}, {2, 3}}, which we write as \ \ \ and \ \ \ . Their normalized 
probabilities are 



Xi 9X3 4 



X\aX_ 



2,3- 



These above two formulas are essentially equivalent to a formula of Kuo [Kuo041 Thms 2.1 
and 2.3]. In the case of In = 6 nodes, we have 



Prl 
Pr( 



1 I 3 I 5 



1 I 3 I 5\ 

2 I 4 I 6 J 



-^3,6(^1,2-^4,5 — -^1,4^2,5) 

-Xi^X 2 ^X 3 ^ + Xi )2 X 34: X 5 q 



and similarly for cyclic permutations of the indices. These cover the five possibilities. 
In the case of 2n = 8 nodes, we have 



P r-r 1 1 3 1 5 1 7 " 

rl l 2 I 8 I 6 I 4 . 



Pn 



1 1 3 I 5 I 7 

2 I 4 I 6 I 8 . 



-(-^1, 2-^7,4 — ^1,4-^7,2) (^3,8-^5,6 — -^5,8-^3,6) 



det 



X\ : 2 


X 1,4 














-^3,6 










-^5,6 




-^7,2 


^7,4 









-^1,2^3,4^5,6-^7,8 + -^1,4X3^X5^X7^ + ^1,6^3,4^5,8-^7,2 + 

-^1,6-^3,8-^5,2-^7,4 + -^1,2-^3,6-^5,8-^7,4 + -X.,4-^3,6-^5, 2-^7,8 — 2X^4X3^X5^X7^ 
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Pv(l\l\l\l)=det 



and similarly for cyclic permutations of the indices. These cover the fourteen possibilities. 

As can be seen from the way we have written these polynomials, for some pairings they 
are expressible as determinants. We plan to discuss such formulas further in [KW08J. 

Of course we could also express the pairing probabilities in terms of the variables X* • = 
Zf J j B /Z WB . The polynomials are exactly the same, although the underlying variables repre- 
sent different quantities. 

1.4 The odd-even pairing projection matrix and explicit formulas 

The explicit computation of the double-dimer pairing probability formulas is quite analogous 
to the computation of the grove partition probability formulas. The role of the L T variables 
for groves is replaced by variables that are indexed by the n! pairings between the odd nodes 
and the even nodes ("odd-even pairings"). If r is such an odd-even pairing, then we define 

X T = X(t) = Y[ ^i,r(i) 
i odd 

and it turns out to be more convenient to work with 

x' T = x\t) = n x i)T(i) 

i odd 

where a cross in a pairing r is a set of four nodes a < b < c < d such that a and c are paired 
with one another and b and d are paired with one another. The U X polynomials" are in fact 
integer linear combinations of the X'^s: 

odd-even pairings r 

As with grove partitions, we construct a projection matrix, but with double-dimer pair- 
ings the matrix projects a vector space with basis vectors indexed by odd-even pairings to 
a vector space whose basis vectors are indexed by planar pairings. Recall that for groves 
the projection matrix has dimensions C n x B n , the n th Catalan number by the n th Bell 
number; for double-dimers the projection matrix p( DD ) has dimensions C n xn\. When n = 4 
the projection matrix T^ 00 ) is 
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We call this matrix p( DD ) the projection matrix from odd-even pairings to planar 
pairings. 

For example, the first row tells us that 



Pr( 



2 14 



yl( 1 I 3 I 5 I 7 ■ 
A \ 2 I 4 I 6 I 8 . 



3 I 5 I 7 ■ 
6 I 2 I 8 . 



-X' 



1 I 3 I 5 I 7 
4 I 8 I 6 I 2 



\r/f 1 I 3 I 5 I 7 \ V ( 1 
- A I 2 I 6 I 8 I 4 J _ A U 
v// 1 | 3 I 5 I 7 \ _ 9 y'( 1 | 3 | 5 ' 
^ I 6 I 4 I 8 I 2 ) I 4 I 6 I 8 

^1,2^3,4^5,6-^7,8 + ^1,2^3,6-^5,8^7,4 + ^1,4^3,6^5,2-^7,8 + -^1,4-^3,8-^5,6-^7,2 

+ -^1,6-^3,4-^5,8-^7,2 ~ *2>X\ 4X3 6 X 5 8 X 7 o + X\ 6 X 3 8 X 5 2 X 7 4. 



, yi( 1 I 3 I 5 I 7 ■ 
~r ^ I 6 I 8 I 2 I 4 , 



As with the matrix V®, we may compute the r th column of the matrix V^ DD ^ via a 
sequence of simple combinatorial transformations on the odd-even pairings. The prototypical 
example of the transformation rule is 



*£4 




4 C 
5 C 



1 1 3 1 5 

4 I 6 I 2 



1 I 3 I 5 1 1 I 3 I 5 1 1 I 3 I 5 
4 I 2 I 6 2 I 6 I 4 6 I 4 I 2 



1 I 3 I 5 

2 I 4 I 6 



1 I 3 I 5 
6 I 2 I 4 



(which is derived in ILemma 3.7ft . Compare with the 16 th column, column \ 
above matrix. More generally, suppose that n > 3 and the odd-even pairing is ^ | ^ 



5 1 : 
2 1 1 

0-2 



in the 
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where the a^'s are odd and the biS are even. Then the transformation rule is 
Z I % I % I rest - H | tl | - | rest + g | « | £ | rest + £ | £ | - | rest 



"3 1 "2 1 "1 1 

- S I g I &r I rest - ^ | ^ | g | rest (Rule 2) 



where in the above "rest" represents £ 4 | • • • \ with unchanged pairings. 

For example, when transforming \ | g 1 1 1 \ we can hold the pair {7, 2} fixed: 

1 I 3 I 5 I 7 ^> 1 I 3 I 5 I 7 _i_ 1 I 3 I 5 I 7 _i_ 1 I 3 I 5 I 7 _ 1 I 3 I 5 I 7 _ 1 I 3 I 5 I 7 
4 I 6 I 8 I 2 — *■ 6 I 4 I 8 I 2 4 I 8 I 6 I 2 8 I 6 I 4 I 2 6 I 8 I 4 I 2 8 I 4 I 6 I 2' 

Of these odd-even pairings, the third and fifth ones are planar, but the others require ad- 
ditional applications of Rule 2. When transforming the first of these terms, if we hold the 
pair {3, 4} fixed, 

1 I 3 I 5 I 7 ^ 1 I 3 I 5 I 7 _i_ 1 I 3 I 5 I 7 _i_ 1 I 3 I 5 I 7 — 1 I 3 I 5 I 7 _ 1 I 3 I 5 I 7 

6 I 4 I 8 I 2 8l4l6l2 "i - el4l2l8 "i - 2l4l8l6 8 I 4 I 2 I e 2 I 4 I 6 I 8 > 

then the resulting odd-even pairings are all planar. The other terms above may be similarly 
transformed into linear combinations of planar pairings, and when they are added up, the 
result is summarized in column \ | g 1 1 1 \ of the projection matrix J>( DD \ 
In ISection 3l we prove 

Theorem 1.4. Any odd-even pairing r may be transformed into a formal linear combination 
of planar pairings by repeated application of \Rule and the resulting linear combination does 
not depend on the choices made when appluina \Rule 2\ so that we may write 

planar pairings a 

For any planar pairing a, these same coefficients vi°^ satisfy the equation 

7WB 

odd-even pairings r 

for bipartite circular planar graphs. 

In ISection 41 we show that the p( DD ) projection matrix of order n is up to signs embedded 
in the V® projection matrix of order 2n: 

Theorem 1.5. 

p (DD) = 

where on the left a and t denote odd- even pairings, and on the right they are interpreted 
as partitions consisting of parts of size 2, and in the sign they are interpreted as maps from 
odd nodes to even nodes, so that o~ x t is a permutation on odd nodes, and (— 1) CT T is its 
signature. 
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1.5 Mult ichor dal SLE connection probabilities 

In ISection 51 we consider the connection probabilities in the scaling limits of the spanning 
tree and double-dimer models, and also of the contour lines in the scaling limit of the discrete 
Gaussian free field with certain boundary conditions. 

Connection probabilities of this sort were first studied by Cardy, who gave (a physics 
derivation for) an explicit formula for the probability (in the scaling limit) of a percolation 
crossing from one segment of the boundary of a domain to another segment of the boundary 
|Car92j . Carleson noticed that Cardy's formula is especially nice when the domain is an 
equilateral triangle, and this was one of the insights that led to Smirnov's proof that the 
percolation interface converges to SLE 6 |Smi01j (see also |CN07j ). (The SLE K process, which 
we do not define here, was introduced by Schramm [SchQOj and describes the scaling limits of 
random curves arising in statistical physics — see [Sch07] .) Arguin and Saint- Aubin [ASA02|, 
§ 3] gave (a physics derivation of) the corresponding crossing probabilities for spins of the 



critical Ising model, shown in Figure 3 As Figure 3 shows, there are two spin interfaces that 



connect the four boundary locations where the boundary conditions change, there are two 
ways in which the interfaces may pair up these boundary points, and the pairing probabilities 
are computed in |ASA02l § 3]. Each spin interface, conditional on the location of the other 
interface, is distributed according to a chordal SLE 3 curve |Smi07j within the domain cut by 
the other interface, so Arguin and Saint- Aubin effectively studied a bichordal version of SLE3. 
Kirchhoff's formula Pr(l|2)/ Pr(12) = can, in view of [LSW04j . be interpreted as giving 
the connection probabilities in a bichordal version of SLE 8 (see § 15. 2p . The corresponding 




Figure 3: The critical 2D Ising model with mixed spin-up and spin-down boundary condi- 
tions, where the spins are shown as black and white hexagons. In this configuration there is 
a connection within the white spins connecting the two white boundary segments. The two 
spin interfaces are shown in bold, and in the scaling limit converge to bichordal SLE 3 . 
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connection probabilities for other values of k were given in [BBK051 § 8.2 and 8.3] (see also 
|Dub06l § 4.1]). 

In a related vein, there has been much study of the physics of multiple-strand networks of 
self-avoiding walks (polymers), loop-erased random walks, and other random paths that are 
either known or believed to be related to SLE K . For example, the "watermelon exponents" 
describe the scaling behavior of the partition function when multiple strands meet at a point 
in the interior or on the boundary of a domain |Sal86[ IDS87[ |Dup92[ |"Dup87| . The boundary 
watermelon exponents for L strands can be extracted from the limiting behavior of the 
connection probabilities of L strands of (multichordal) SLE connecting 2L boundary points 
when L of the boundary points are close to each other and the other L boundary points are 
close to each other, though of course these connection probabilities contain more information 
than is summarized in the exponents. This approach does not yield the bulk watermelon 
exponents, since to obtain the bulk exponents one would need the endpoints of the strands 
to be in the interior of the domain. The scaling exponents associated with more complicated 
multiple-strand networks (for 2 < k < 8) have also been derived |Dup87 Dup89 , IOB88] . See 
|Dup06| for up-to-date lecture notes on the physics of networks of polymers and other types 
of strands and their relation to SLE. 

It is only natural to consider the connection probabilities of multichordal SLE K , which 
is defined in |KL07j . and has the property that, conditional on one curve, the remaining 
curves are distributed according to multichordal SLE K with fewer curves within the domain 
cut by the conditioned curve. Kozdron and Lawler [KLQ7] studied multichordal SLE K for 
a fixed connection topology (when k < 4). Cardy exhibited several discrete models with 
multiple interfaces that arise naturally in the context of conformal field theory, for which 
the interfaces have the same joint distribution (whose scaling limit is thought to be multi- 
chordal SLE) conditional on a particular connection topology, but for which the connection 
topologies have different probabilities [Car 07j. (Thus to discuss connection topology proba- 
bilities of multichordal SLE, one must specify a discrete model whose scaling limit is being 
taken.) Comparing two such models, the connection topology probabilities in one model 
have different weights compared to another such model, but these weights do not depend 
on the domain or on the location of the boundary points where the boundary conditions 
change, so that the set of all connection topology probabilities in one such model determines 
the connection topology probabilities in another such model. Dubedat [Dub06] analyzed the 
general multichordal SLE K connection probabilities (and included special discussion of the 
cases k = 2,6,8), but there still remain open problems regarding these connection probabil- 
ities when the number of curves is large. 

The scaling limit calculations in ISection 51 yield explicit formulas for the multichordal 
SLE K connection probabilities in the cases k = 2, 4, 8 (for any number of curves, for any 
connection topology). There is a scaling limit relation |LSW04j between branches of uniform 
spanning trees on periodic planar graphs and SLE 2 . Essentially, the scaling limit, as the 
lattice spacing tends to zero, of a branch of the uniform spanning tree on a bounded domain, 
tends to a random simple curve which is equal in law to SLE2. Similarly, the curve which 
winds between the uniform spanning tree and its dual spanning tree (|Figure 9|) converges 
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in the scaling limit to an SLEg [LSW04]. The double-dimer paths are thought to have a 
scaling limit that is given by SLE4, but this has not been proved. However, the contour 
lines in the discrete Gaussian free field with certain boundary conditions have been proved 
to converge to SLE 4 [SS06] in the scaling limit, and using the results of [SS06J we prove in 
ITheorem 5.1l that the probability distribution of the pairings of these contour lines coincides 
with the pairing distribution for the double-dimer model. 



2 Grove partitions 



2.1 The meander matrix 

Associated to a planar partition <j on n nodes is a planar chord diagram consisting of n 



disjoint chords, winding between the components and the dual components of a, see |Figure 4 
for an example. The chords have in total 2n endpoints, one between each node and adjacent 
dual node. A planar chord diagram determines a and vice versa. 





1 


3 


5 


7 1 9 1 11 


13 


15 


2 


16 


10 


6 1 8 1 12 


4 


11 



planar chord diagram 



There is a natural bilinear form on the space of formal linear combinations of chord 
diagrams that is defined as follows. Given two chord diagrams G\ and C2, we can draw them 
on a sphere, one in the upper hemisphere and one in the lower hemisphere, so their common 
boundary consists of 2n points around the equator. In the resulting figure the chords join up 
to form simple closed loops, each crossing the equator (at least twice) and such that the set 
of all loops crosses the equator a total of 2n times. Such an object is a meander of order 
n. The bilinear form (Ci,C2)q is defined by 



(C u C 2 ) 9 = q 



# loops in meander formed from Ci and C2 



(2) 



This definition can then be extended linearly to formal linear combinations of chord diagrams. 

The matrix A4 q which has rows and columns indexed by chord diagrams and has en- 
tries (■M q )c 1 ,c 2 — {Ci,C2) q is known as the Gram matrix of the Temperley-Lieb algebra 
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(also called the meander matrix). Jones [Jon83j determined the values of q for which A4 q 
is positive definite and positive semidefinite, Ko and Smolinsky [KS91] determined when 
A4 q is nonsingular, and Di Francesco, Golinelli, and Guitter [DFGG97] gave an explicit di- 
agonalization of M. q . For our results on groves we shall use the fact that \im q ^ A4 q /q is 
nonsingular |DFGG97l Eqn (5.18)], and for our results on the double-dimer model we shall 
use the fact that M. 2 is nonsingular [DFGG971 Eqn (5.6)]. 

2.2 The "L" polynomials 

In this subsection we show that Pr(cr) is a polynomial in the L^j for each planar partition cr. 

For a planar partition a let C a be its associated chord diagram. We define a bilinear 
form (,) t on the vector space generated by planar partitions by 

. . (C ai C T )a I 1, if meander formed from C a and C T has one loop, . . 

(cr, r) t = hm - = < (3) 

<?^° q 10, otherwise. 

We let M. t = lim g ^ M. q /q denote the Gram matrix of (,) t , which is nonsingular. 

For a partition cr on {1, . . . , n} and graph Q with nodes {1, . . . , n}, let Q a be the graph 
Q with nodes identified according to cr, that is Q a = Q / ~ where v ~ v' if and only if v and 
v' are in the same component of a. 

If cr and r are planar partitions and Q is a circular planar graph, it is not hard to see 
that 

{1, if a grove of type r in Q gives a spanning tree in Q a 
when the vertices are identified according to a, (4) 
0, otherwise. 

Indeed, given a grove of Q of type r, and a "grove" of type a living within the outer face of Q, 
if their union is a spanning tree, the path winding around the outside of this spanning tree 
is a meander. We can contract to points the components of the "grove" of type a, without 
changing the topology of the meander. Another way to say this is that (cr, r) t = 1 if and 
only if the number of parts of a and r add up to n + 1 and the transitive closure of the 
relation defined by the union of a and r has a single equivalence class. 

Let Z(t) be the weighted sum of groves of type r on Q and Zg^r) the corresponding 
weighted sum on Q a . In particular Zg CT (tree) is the weighted sum of spanning trees of Q a . 

Lemma 2.1. For a planar partition a and circular planar graph Q, 

T 

where the sum is over all planar partitions. 
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Proof. If (a, r)t = 1 then every grove of Q of type r, when we identify vertices in a, becomes 
a spanning tree of Q a . Conversely, any spanning tree of Q a arises from a unique grove of Q, 
and this grove has partition r satisfying (a, r) t = 1. □ 

Note also that Zg a (uncrossing) = Z (uncrossing), since an uncrossing of Q is the same as 
an uncrossing of Q a . 

Lemma 2.2. For any graph Q (not necessarily planar) and partition o (not necessarily 
planar), Zg a (tree) jZ (uncrossing) is a polynomial in the L^j. 

Proof. Let S be the matrix whose rows index the equivalence classes of the relation a and 
whose columns index the nodes 1, . . . , n, and whose i, j-entry is S^j = 1 if node j is in class i, 
and Sij = otherwise. 



The response matrix A(Q a ) for Q a (see Appendix A ) is obtained from the response ma- 
trix A for Q simply as 

A(&) = SAS T . (5) 

That is, putting potential Vi on vertex i of Q a is the same as putting potential Vi on each 
node of Q in the equivalence class i. The resulting current out of node i of Q a is the sum of 
the currents out of nodes of Q in the i th equivalence class. In particular, entries in A(Q a ) are 
sums of entries in A. From [Lemma A.TI we have 

det A(Q a ) = Z Qa (tree) jZ Qa (uncrossing), 

where A{Q a ) is A(Q a ) with one row and column removed. The left-hand side is a polynomial 
in the L it j, by (jSJ), and since Zg a (uncrossing) = Z(uncrossing), this concludes the proof. □ 

Using [Lemma 2.11 we can invert matrix M. t (since it is nonsingular) to write Z(r) as a 
rational linear combination of the Zg a (tree) as a varies. Dividing both sides by Z(uncrossing) 
we see that Pr(r) is a polynomial in the Ljj with rational coefficients. In the next subsection 
we prove that (even though .M^ 1 has noninteger entries) the coefficients of these polynomials 
are actually integers. 

2.3 Integrality of the coefficients 

We can extend the bilinear form (,) t to work on any pair of partitions (not necessarily planar), 



simply by taking the characterization in Equation 4 and dropping the requirement that the 



partitions be planar. Define an "extended meander matrix" £ t with rows indexed by planar 
partitions and columns indexed by all partitions as above, so that {S t ) UtT = (cr, r) t . The 
matrix S t has dimensions C n x B n , where C n is the n th Catalan number, and B n is the n th 
Bell number (the number of partitions on n items). The matrix M. t is the submatrix of S t 
containing only the columns for planar partitions. 

Note that with this extended definition of (,) t , ILemma 2.11 holds for general, nonplanar 
graphs, provided we sum over all (not necessarily planar) partitions. We let G be the 
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column vector of "glue variables" , whose entries are G a = Zg a (tree) for a a planar partition 
of 1, . . . , n. Let Z be the column vector of partition variables, whose entries are Z T = Z(t) 
where r runs over all partitions. The extension of lLemma 2.11 to the nonplanar setting gives 

G = StZ. (6) 

For any not necessarily planar graph Q on n nodes there is an electrically equivalent 
complete graph JC on n vertices (in which every vertex is a node): it is the graph whose 
edge {i,j} has conductance L it j(Q, N). (The graphs Q and K are electrically equivalent in 
the sense that, when the same voltages are applied to the nodes of Q and /C, the current 
responses will be the same, i.e., they have the same Dirichlet-to- Neumann matrix.) Note 
that each Z^ir) is trivially a polynomial in the L it j, since each crossing of K, has weight which 
is a monomial in the Ljj. In fact, we can write this polynomial explicitly: for each part A of 
r, we count the weighted sum of spanning trees of the complete graph on the vertices in A, 
and then take the product over the different parts of r. Recalling our definition of L T in 



Equation 1, we have that Zjc(t) = L T . For example 

Zk:(1|23|456) = 1 ■ ^2,3 ■ (£4,5-^5,6 + £4,5-^4,6 + -^4,6-^5,6) = -^1|23|456- 

We let L be the column vector of partition variables for /C: L T = Zjc(t) = L T . 



By the preceding equation, for the graph K,, Equation 6 specializes to G = £ t L. By 



ILemma 2.2\ for any graph, G is determined by the L's, so we have 

G = £ t Z = S t L 

for any graph. For planar graphs, the entries of Z corresponding to nonplanar partitions are 
0, so £ t Z = M.tZ, and hence M. t Z = £ t L. Since Ai t is invertible, 

Z = M~ x E t L. 

We define 

We shall see how to compute directly, i.e. without inverting A4 t . The direct computation 
of V® involves only integer operations, from which it will follow that the "L polynomials" 
have integer coefficients. 

What is the matrix pW? For each planar partition a, the row a tells usZ ff = ^ t V®L t . 
For each partition r, the r th column gives us a linear combination of planar partitions 
XL r&l-o that is equivalent to r in the sense that for any planar partition p, 
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(We shall see that <^p, Pi+ery = (p, r) t for nonplanar partitions p too.) 

Let us say that two linear combinations of partitions on n items a T r and Y2 T Pt t 

are equivalent (=) if for any (possibly nonplanar) partition p on n items J2 T a r(P^ T )t — 
J2 T Pr(p,T) t . For example, 

Lemma 2.3. 1|234 + 2|134 + 3|124 + 4|123 = 12|34 + 13[24 + 14|23 

Proof. For any partition p with three parts, (p, LHS) t = 2 = (p, RHS) t : By symmetry 
considerations, we need only consider one such parition, say 12|3|4, and 

(12|3|4, LHS) t = 1 + 1 + + = + 1 + 1 = (12|3|4, RHS) t . 

For partitions p with other numbers of parts, (p, LHS) t = = (p, RHS) t , since from 
Equationl] (p, r) t = whenever #{parts of p} + #{parts of r} ^ #items + 1. □ 

As we shall see, this lemma, together with the following two lemmas, which show how 
to adjoin new parts and new items to the partitions of the left- and right-hand sides of an 

equivalence (=), will allow us to write any partition as an equivalent (=) sum of planar 
partitions. 

Lemma 2.4. Suppose n > 2, r is a partition of 1, . . . , n — 1, and r = a a a. Then 

r\n = a a a\n. 

17 

Proof. If {n} is a part of p then (p, r|n) t = = (p, RHS) t . Otherwise (p, r|n) t = (p\n, r) t = 

If t is a partition of 1, . . . , n — 1 and j G {1, . . . , n — 1}, we can insert a new item n into 
the part of r that contains item j. We refer to the resulting partition on 1, . . . , n as "r with 
n inserted into j's part." 

Lemma 2.5. Suppose n > 2, r is a partition of 1, . . . , n — 1, j G {1, . . . , n — 1} ; and 
r = J2 a Oi a a. Then 

[t with n inserted into fs part] = ^^a CT [cr with n inserted into j's part]. 

Proof. If j and n are in the same part of partition p, as well as being in the same part 
of partition it, then by (TJJ (p, 7r) t = 0. Thus (p, [r with n inserted into j's part]) t = = 
(p, RHS) t . If j and n are in separate parts of p, then let p' denote the partition obtained 
from p by merging the two parts containing j and n, and then deleting n. By (J3J) we have 

(p, [r with n inserted into j's part]) t = (p, r) t = a a {p\ cr) t = (p, RHS) t . □ 
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Lemmas 12.31 12.41 and 12.51 imply that the left-hand-side and right-hand-side of transfor- 
mation Rule 1 are equivalent (=). 

Theorem 2.6. For any partition t, there is an equivalent (=) integer linear combination of 
planar partitions J2a a ^ a (i- e -> where each a a G Zj. 

Proof. We prove this by induction on the number of items in the partition. The theorem is 
true for n < 3 since each such partition r is already planar (and with lLemma 2.3[ we see that 
it is also true for n = 4). Suppose r contains more items. If {n} is a part of r, then we may 
use the induction hypothesis together with lLemma 2 .41 to find the desired linear combination 
of planar partitions. Otherwise, item n is in the same part as some other item j (if there 
is more than one choice of j, it does not matter which one we pick). By the induction 

hypothesis, we may write r \ n = an integer linear combination of planar partitions, and by 
ILemma 2.51 we may write r as an integer linear combination of "almost planar" partitions, 
by which we mean partitions that would be planar if the item n were deleted from them. 

Next we use Lemmas l2.3[l2.4l[2.5l to express an almost-planar partition /i as an equivalent 
integer linear combination of planar partitions. We shall use induction on the number k of 
parts of that cross the chord from j to n. There is nothing to show if A; = 0, and 
otherwise we consider the part S of // closest to j that crosses the chord from j to n. 
We let a = {i G S : i < j} and c = {i G 5* : i > j}, both of which are nonempty, 
b = the part of u. containing j, and d = {n}. Let oq G a and Cq G c. From [Lemma 2.3"! upon 
relabeling 1 — > a , 2 — > j, 3 — > Co, and 4 — > n, we get 

a ,c \j,n = a \j,c ,n + j\a ,c ,n + c \a ,j,n + n\a ,j,c - a ,j\c ,n- a ,n\j, cq. 

Using [Lemma 2 .51 we may insert the rest of a into the parts containing ao, the rest of c into 
the parts containing cq, and the rest of b into the parts containing j, to obtain 

a U c\b U d = a\b U c U d + b\a U c U d + c\a U b U d + d\a U b U c - a U b\c U d - a U d\b U c. 

By further application of Lemmas 12.41 and 12.51 we may adjoin each of the remaining parts of 
H to the left-hand-side (thereby obtaining a) and to each partition on the right-hand-side. 
Of the resulting partitions on the right-hand-side, the fourth one is planar, and the other 
partitions are almost planar with k — 1 parts crossing the chords from n to the rest of n's 
part. □ 

For a partition r, let ^ CT a a a be any linear combination of planar partitions equivalent 
(=) to t. Now J2a(^^ T ~ a cr)< J lies in the null-space of A4 t , but A4 t is nonsingular, so 
Va)r = tt<r for each a. In particular, the linear combination promised by ITheorem 2.61 is 
unique, and gives the r th column of V® ■ This linear combination was obtained by repeated 
application of Rule 1, which completes the proof of ITheorem 1.21 which in turn implies the 
second part of ITheorem 1.11 as a corollary. 

We remark that the entries of the matrix are all or ±1 for n < 7 nodes, but that 
when n > 8 other integers appear. 
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2.4 The U R" polynomials 

In this subsection we prove the first part of ITheorem 1.1 1 i.e., we show that Pr(cr) is a 
certain polynomial in the Rij's for each planar partition <j. We start with some elementary 
characterizations of the Rij and Ly variables in terms of groves. 

Proposition 2.7. The resistance between i and j is 

A:i<=A,j£A c 

Proof. If nodes i and j were the only two nodes, then we know from Kirchhoff 's formula (the 
case n = 2) that Ri j is the ratio of the probability that i and j are in different components 
of a 2-component random grove, to the probability of a 1-component grove. When there 
are other nodes, these 2-component groves necessarily take the form A\A C where i G A and 
j E A c . (See the top half of |Figure~6j ) □ 



The following proposition may be derived from Lemma 4.1 of [CIM98j . but for the reader's 
convenience we provide a short proof. 

Proposition 2.8. For i ^ j , 

Lij = Pr(z, j | rest singletons), 

i.e., Lij = Pr(er) where a is the partition in which every part is a singleton except for the 
part {i,j}. 

Proof. Recall from the construction in lTheorem 2.61 that whenever a partition r is expressed 
as an equivalent sum of planar partitions, each of the partitions has the same number of parts 
as r. Since o has n — 1 parts, and any partition r with n — 1 parts is already planar, the 
cr th row of is nonzero only in the a th column, so Pr(cr) = L a = Lij. □ 



Any circular planar graph Q has a dual circular planar graph Q* , as shown in |Figure 5 
The dual Q* has an inner vertex for every bounded face of Q, and a node numbered i between 
consecutive nodes i, i + 1 mod n of Q. For each edge of Q there is a dual edge of Q*, whose 
conductance is, by definition, the reciprocal of the conductance on the corresponding edge 
of Q. For each grove of Q there is a dual grove of Q* formed from the duals of the edges of 



Q not contained in the grove (see Figure 5). A grove in Q has weight equal to the weight of 
its dual grove in Q*, times the product of the conductances of all edges in Q. 

Proposition 2.9. For the dual graph we have 
and 

"I, EE ( 8 ) 

\<i'<j'<n 
chord i'.j' crosses dual chord i,j 
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Figure 5: Shown here is a circular planar graph Q with four nodes (shown in black) and 
one inner vertex, with edges shown as solid lines. The dual circular planar graph Q* has 
four nodes (shown in white) and two inner vertices, with edges shown as dashed lines. Also 
shown is a grove of Q of type 1 1 234 (the edges of the grove are shown in bold) and its dual 
grove of Q* (with edges also shown in bold), which has type 14|2|3. 



Proof. By Proposition 2.7 



7?* 

H i,3 



^A-i&A jeA c P r (^l^ c )- But dual nodes i and j are in 
opposite parts of A\A C if and only if the partition dual to A\A C is i', j'jrest singleton for some 
chord i',f crossing dual chord i,j (see |Figure 6\j; applying proposition 2.8|then yields (JSJ). 



Expanding the right-hand-side of ([7j) using (jSJ) yields the left-hand-side of ([7]). 



□ 




Figure 6: 4 can be expressed as a sum of 2-part partitions in which 1 and 4 are in different 
parts. For the dual resistance R\ 4 in the dual graph, the sum over dual partitions becomes a 
sum over partitions consisting of a doubleton part and rest singletons, where the doubleton 
part separates dual nodes 1 and 4. 
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Proof of first part of \Theorem Since the dual of a partition a is a partition a* on the 
dual graph, and each grove has the same weight as its dual grove, times a constant, we have 

Pr(o-) Pr(a*) 



Pr(tree) Pr (dual uncrossing) ' 

which by the second part of ITheorem 1.11 is an integer-coefficient polynomial in the L*j's, 
which by Proposition 2.9 is an integer-coefficient polynomial in the -Rjj/2's. □ 



3 Double-dimer pairings 



Recall that for our double-dimer results we assume the circular planar graph Q is bipartite, 
so we may color its vertices black and white so that each edge connects a black vertex to 
a white vertex. It is convenient to assume that the nodes on the outer face alternate in 
color, so that the odd-numbered nodes are black and the even-numbered nodes are white. 
If the graph Q does not satisfy this property, we can extend Q by adjoining an extra vertex 



and edge with weight 1 for each node that has the wrong color (refer to Figure 7 ), and the 
double-dimer configurations of the extended graph are in one-to-one weight-preserving and 
connection-preserving correspondence with the double-dimer configurations of the original 




Figure 7: Double-dimer configuration on a rectangular grid with 8 terminal nodes. In this 
configuration the pairing of the nodes is {{1, 8}, {3, 4}, {5, 2}, {7, 6}}, which we write as 
g | 4 | 2 | 6 • The double-dimer configurations on the graph on the left are in one-to-one 
correspondence with the double-dimer configurations of the extended graph on the right, for 
which the odd terminals are colored black and the even terminals are colored white. 
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graph. Furthermore, the quantities Z , Z as well as the variables Xij are the same on 
this new graph as they were on Q. Henceforth we assume without loss of generality that the 
node colors of Q alternate between black and white. 

3 . 1 Polynomiality 

Since the nodes alternate in color, Z BSN is the weighted sum of dimer covers of Q. Recall that 
Z DD is the weighted sum of double-dimer covers of (£?, N) where each of the terminal nodes 
is included in only one of the dimer coverings. Let S be a balanced subset of nodes, that is, 
a subset containing an equal number k of white and black nodes. Let Z D (S) = Z D (Q \ S) 
be the weighted sum of dimer covers of Q \ S, and Z D = Z D (0) = Z BW . The superposition 
of a dimer cover of Q \ S and a dimer cover of Q \ S c , where S c = N \ S, is a double-dimer 
cover of (Q, N). In fact, we have 

Lemma 3.1. Z D (S)Z D (S C ) is a sum of double-dimer configurations for all connection topolo- 
gies tt for which tt connects no element of S to an element of S c . That is, 



where Ms >7T is or 1 according to whether tt connects nodes in S to S c or not. 

As a special case, when S = we have Z D (0)Z D (N) = Z DD , which is closely related 
to Ciucu's graph factorization theorem [Ciu97t Thm. 1.2]. Ciucu showed how to enumerate 
dimer coverings in a bipartite graph with bilateral symmetry. Dimer coverings of the whole 
graph correspond to double-dimer configurations of one half of the graph, with vertices on 
the symmetry axis corresponding to nodes, except that there is an extra factor of 2 in weight 
for each path connecting a pair of nodes, unless the path consists of a single edge. So the 
graph factorization theorem contains a power of 2 for each pair of nodes, and a weight of 1/2 
for each edge on the symmetry axis, neither of which appear in ILemma 3.11 

Our remaining double-dimer results only make sense if there are double-dimer configu- 
rations of Q, which implies that there are dimer configurations of Q, so that Q has equal 
numbers of black and white vertices. After proving this lemma, we shall henceforth make 
this assumption. 

Proof of \Lemma 3.1\ Each double-dimer path connecting a pair of nodes separates an even 
number of nodes on its two sides, and since the white and black nodes alternate around the 
boundary, each double-dimer path must go from a node to a node of the opposite color. 

Consider the double-dimer cover of (Q, N) formed by the superposition of a dimer cover 
of Q \ S and a dimer cover of Q \ S c . On the double-dimer path starting from a white node 
of S, every dimer from a white vertex to a black vertex is from the second dimer cover and 
every dimer from a black vertex to a white vertex is from the first cover. So such a path 
necessarily ends at a black vertex in S. Similarly a path from a black vertex in S necessarily 
ends at a white vertex in 5*. 




(9) 



TT 
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Conversely, any double-dimer configuration with no connections from S to S c can be 
decomposed into a dimer cover of Q \ S and a dimer cover of Q \ S c . There are 2 e choices of 
such a decomposition, 2 choices for each of the £ closed loops. □ 

Before proceeding with the proofs of Theorems 11.41 and 11.31 we recall some basic facts 
about Kasteleyn matrices. Kasteleyn matrices may be used to enumerate dimer coverings of 
any planar graph |Kas67j . but we review here just the bipartite case. A Kasteleyn matrix 
of an edge-weighted bipartite planar graph (with a given embedding in the plane) is defined 
to be a signed adjacency matrix, K = {K w ^) with rows indexed by the white vertices and 
columns indexed by the black vertices, satisfying the following: K w b is ± the weight on edge 
wb and if there is no edge. The signs of the edges are chosen so that around each face there 
are an odd number of — signs if the face has mod 4 edges and an even number of — signs 
if the face has 2 mod 4 edges. Kasteleyn |Kas67j showed that every bipartite plane graph 
with an even number of vertices has a Kasteleyn matrix, and if there are equal numbers of 
black and white vertices then | deti^| is the weighted sum of all dimer coverings, where the 
weight of a dimer covering is the product of its edge weights. 

From Kasteleyn theory [Kas67j [Ken97j it is straightforward to compute Z D (S), though 
some work is needed to get all the signs right. But the signs for Z D (S)Z D (S C ) are simpler, 
and this product is all we need anyway. Recall that Xij = Z D ({i, j})/ Z D . 

Lemma 3.2. Let S be a balanced subset of {1, . . . ,2n}. Then 

Z D (S)Z°(S°)/(Z°T = det [(l i>jes + l lMS ) x i-ir-^X^^X 1 ■ ( 10 ) 

Proof. For convenience we adjoin to the graph Q 2n edges along the outer face connecting 
adjacent terminal nodes, and give these edges weight (or weight e and then take the 
limit e I 0). Given a Kasteleyn matrix of a graph, the signs of edges incident to a vertex 
may be reversed, and each face will still have a correct number of minus signs. For each 
i — 1, . . . , 2n— 1, in order, if the edge from node i to node i + 1 has a minus sign, let us reverse 
the signs of all edges incident to node i + 1. Doing this ensures that for 1 < i < 2n — 1, the 
sign of the edge from node i to node i + 1 is positive. The sign of the remaining adjoined 
edge, from node 2n to 1, will necessarily be — (— l) n for the outer face to have a correct 
number of minus signs. 

Let (w\, bi), . . . , (wk, bk) be any noncrossing pairing of the nodes of S, where Wi, . . . , io& 
are the white nodes of S and b\, . . . ,bk are the black nodes of S. Let us adjoin edges of 
weight W to the outer face connecting W{ to bi for 1 < i < k. To retain the Kasteleyn 
sign condition, the sign of a new edge connecting black node b and white node w will be 
(— l)(l b ~ w 'l~ 1 )/ 2 . Let K\y be the Kasteleyn matrix of the resulting graph, with the rows and 
columns ordered so that W\, . . . , are the first k rows and b\, . . . , b}. are the first k columns, 
and let K = Kq be the corresponding Kasteleyn matrix when W = 0. Recall that [x Q ]p(x) 
denotes the coefficient of x a in the polynomial p(x). Then Z D (S) = ±[W k ] det Kw and 
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Z D = ± det Kq. But det Ky/ enumerates the weighted matchings of the enlarged graph with 
all the same sign, so 

Z^S) [W k ] tetK w ^tobl det^ Eti Mbi x 6=6l) ... i6fe 

Z D " [W°]detAV " detK 1 ; " 4 M»J«=«i.-..,«* 

(11) 

where lf\g denotes the sub matrix of formed from deleting the rows and columns from 5*, 
and the last equality is Jacobi's determinant identity. The special case S = {b, w} yields 

X b , w = Z D ({b,w})/Z D = (-l)(l^l-D/2^-i_ 

From this equation and (ITT1) we get (flOl) . except possibly for a global sign change. 

Next we compare the sign of Z D (S)Z D (S C )/(Z D ) 2 that we get from (fTTj) to the sign 
in PH]) . In the event that 5 = {1, . . . , 2/c}, we can take b±, . . . , bk = 1,3, ... ,2k — 1 and 
W\, . . . ,Wk = 2,4,..., 2k, and for S c we can take bi, . . . , b n -k = 2k + 1, 2k + 3, . . . , 2n — 1 
and w x , ... , w n _ k = 2k + 2, 2k + 4, . . . , 2n, so we see that the sign for Z D (S)Z D (S C ) / \Z D f 
from (TTTj) agrees with (flOj) . Now suppose that for some S the signs from (fTT!) agree with (JTUJ) , 
and we replace one of the nodes s 6 5 with s + 2 ^ S to get S". If s + 1 ^ S, the power 
of —1 in (flTi) changes by 1. If the s + 1 £ 5, we may assume s was paired with s + 1, 
and we see that the power of —1 in (1111) does not change. Thus the power of —1 from (1111) 
is opposite for Z D (S')Z D (S' C )/(Z D ) 2 compared to Z D (S)Z D (S C )/(Z D ) 2 . But the product of 
determinants from (fill) is ± the determinant from (TlO|) . and the choice of sign is opposite 
for (S, S c ) and (S", 5" c ). Thus ( TTOl) has the correct global sign for 5", and hence by induction 
for any balanced set of nodes. □ 

By ILemma 3.21 the quantities Z D (S)Z D (S C ) /(Z D ) 2 are homogeneous polynomials of de- 
gree n (in fact determinants) in the X^j, so the matrix M from ([9]) maps the quantities 
Pr(7r)Z DD /(Z D ) 2 to homogeneous polynomials of degree n in the Xij. We need to show that 
the matrix M has full rank, that is, the rank of M is C n , the n th Catalan number, which is 
the number of different planar partitions 7r. 



Lemma 3.3. M T M = M.2, the meander matrix M q from (T5]) \DFGG97l evaluated at q = 2. 



Proof. Given two planar matchings tt, a, the integer S a M T M5 n is the number of subsets S 
with the property that neither tt nor a connects a node in S to a node in S c . 

We draw it and o on the sphere, with tt in the upper hemisphere and a in the lower 
hemisphere, so that their union is a meander crossing the equator 2n times. Suppose this 
meander has k components. For each component, the nodes alternate black and white, and 
so the component has the same number of black nodes as white nodes. For each component 
we choose to put all of its nodes in S or all of its nodes in S c . We can construct 2 k possible 
sets S in this way, and these are exactly the sets S for which neither tt nor a connects a 
node in S to a node in S c . 

Thus an entry in M T M which corresponds to a meander with k components is 2 k . □ 
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By [DFGG97] Eqn (5.6)] the determinant of the meander matrix at q = 2 is 

n 

i=i 

where certain integers. In particular it is nonzero, so M has full rank and we 

can solve ® for Pr(7r)Z DD /(Z D ) 2 in terms of the quantities Z D (S)Z D (S C )/(Z D ) 2 which by 
ILemma 3.21 are homogeneous polynomials in the X^j. This effectively proves ITheorem 1.3} 
except for the part about the coefficients being integers. 



3.2 Integrality of the coefficients 

We start by collecting the vectors and matrices we need. Extend the Gram matrix M T M 
when q = 2 to a C n x n\ matrix £2 whose rows are indexed by planar pairings and whose 
columns are indexed by all (not necessarily planar) pairings connecting odd nodes to even 
nodes. We have the following matrices and vectors. 

M =matrix from balanced subsets S C {1, . . . , 2n} to planar pairings n, from 
M.2 =M T M = the C n x C n meander matrix (with q — 2) for planar pairings 
£2 =C n x n\ extended meander matrix (with q = 2) for odd-even pairings 
P =vector of normalized pairing probabilities P w = Pr(n) = Pr(vr)Z D ({l, . . . ,2n})/Z D 

indexed by planar pairings 7r 
D =MP = vector of products D s = Z D (S)Z D (S C )/(Z D ) 2 

indexed by balanced subsets S C {1, . . . , 2n} 
X =vector of X-monomials indexed by odd-even pairings; X p = Yl{ij}ep-^hi 
X' =vector of X-monomials as above with sign (-1)* crosses of p 

Recall that we defined a cross of a pairing p to be a set of two parts {a, c} and {b, d} of p 
such that a < b < c < d. We define the parity of an odd-even pairing p = ^ | ^ | • • • 1 2 " _1 to 
be the parity of the permutation (wi/2)(w 2 /2) . . . (w n /2). We shall use the following fact 

Lemma 3.4. For odd-even pairings p, 

^_^parity of p J~J ^_^(|i-j|-l)/2 _ crosses of p j-^) 

(*,i)ep 

Proof. When p = \ \ \ \ ■ ■ ■ both sides of the above equation are +1. Now suppose we do a 
transposition, swapping the locations of w and w + 2. (Such swaps are enough to connect 
the set of odd-even pairings.) In the event that one of w or w + 2 was paired with w + 1, 
such a swap will not change the sign of the left-hand-side (it changes the sign of exactly 
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one term in the product and also the parity of p), nor does it change the number of chords 
that cross. Otherwise the sign of the left- hand- side does change (exactly two terms of the 
product change sign, and the parity of p also changes). Also, the number of chords crossing 
w + l's chord changes by or ±2, and the chords containing w and w + 2 now cross if they 
didn't before, and vice versa. □ 

Lemma 3.5. M T D = £ 2 X' 

Proof. Recall from ILemma 3.21 that 

Ds = ZD( g? C) = det [(l tJes + 1,,^) x (-l)''--'- 1 '^,^:™ • 
When we expand the determinant, we get 

D S = J2 (_l)Parityofp "Q (_l)([^i[-l)/2jjf y 

odd-even pairings p (*j)£p 
p does not bridge .S' to S c v / 

= X' p by G2J 

Let 7r be a planar pairing. Upon summing over sets S that are not bridged by it we get 



odd-even pairings p, 
7r, p do not bridge S, S c 

E E K 

Id-even S:ir,p do not 
irings p bridge S 

^ ^ 2# cycles in p U 7r y^l 



S'C{l,...,2n} odd-even pairings p, 

7r does not bridge S to S c S: it, p do not bridge S, S' 



odd-even 
pairings p 

The left-hand side is the 7i~ th row of M T D, and the right-hand side is the n th row of EzX' . □ 
Theorem 3.6. A^ 2 ^ = ^X' 

Proof. Since MP = D, we have A4 2 P = M T MP = M T D = £ 2 X'. □ 
Since A^2 is invertible, we may define 

p( DD ) = M 2 %. 

Since P = V^ DD ^X', we can interpret p( DD ) as the matrix of coefficientsjjf the X' polynomials: 
for a given planar pairing cr, the a th row of p( DD ) gives the polynomial Pr(cr). The r th column 



Boundary partitions in trees and dimers 



Kenyon <fc Wilson 



27 



of V {DD) gives a linear combination of planar pairings that is in a sense equivalent under (,)2 
to r: for any planar pairing p, 



\ c I 2 <r 



(We shall see (p, ^2 a V^^ aj = (p, r) 2 for nonplanar odd-even pairings p too.) 

We say that two linear combinations of odd-even pairings a T r and /3 T r are equiv- 
alent (=) if for any (not necessarily planar) odd-even pairing p, we have ^2 T ct T (p 1 T} q = 
Y2 T Pt(P) T )q- Following the approach we used in lSection 2.31 for groves, we will show how to 

2 

transform an odd-even pairing r into an equivalent (=) linear combination of planar pairings, 
and thereby compute p( DD ) using only integer operations, i.e. without inverting M.2- 
The following key lemma is analogous to ILemma 2.31 



T nmm o Q 7 1 I 3 I 5 i 1 I 3 I 5 i 1 I 3 I 5 1 I 3 I 5 I 1 I 3 I 5 I 1 I 3 I 5 
Ijemilld O.i . 2 I 4 I 6 4 I 6 I 2 6 I 2 I 4 — 4 I 2 I 6 ' 2 I 6 I 4 i" 6 I 4 I 2 



Proof. For any odd-even pairing p, we have (p, LHS)2 = 12 = (p, RHS)2- For example, if 
p = then (p,LHS) 2 = 2 2 + 2 2 + 2 2 = 12, while (p,RHS) 2 = 2 3 + 2 1 + 2 1 = 12. 

For the other odd-even pairings p on {1,...,6}, it is similarly straightforward to verify 
(p, LHS) 2 = (p, RHS) 2 because 2 1 + 2 1 + 2 3 = 2 2 + 2 2 + 2 2 . □ 

ILemma 3.71 is a special case of ILemma 3.81 below. Though we do not need ILemma 3.8f s 
extra generality for our results, some readers may find its proof more instructive. 

Lemma 3.8. If n is a positive integer, then 

E (-i)v^o 

odd-even pairings a on 2n items 

if and only if q is an integer satisfying < q < n. 

Proof. If cr is a permutation on {1,2, ...,n}, we can interpret a as the odd-even pairing 
2ax 1 2<J2 1 ' ' ' 2™ n _1 an d vice versa. Note that if p and a two odd-even pairings, then the cycles 
in the union of p and a are the cycles of p _1 cr when p and a are interpreted as permutations, 

i.e., (p,a) q = g# cycles in p-! CT _ Note algo that (_!y = cycles in p^a Thug 

p, J2 (-i)^/ =(-i) n (-irEc(n,fc)(-# 

odd-even pairings cr on 2n items / k=l 

where the c(n, k) are the unsigned Stirling numbers of the first kind, which count the number 
of permutations on n letters that have k cycles. It is well known (see e.g. |Sta86| Proposi- 
tion 1.3.4]) that J2k c (™> k ) xk =x(x + l)(x + 2)---(x + n-i). □ 
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The next lemma is analogous to Lemmas 12.41 and 12.51 



Lemma 3.9. Suppose n >2, r is an odd-even pairing of 1,2, ... , 2{n — 1), and r = a o-°~- 
Then 



i2n-l X n ~|2n-l 

' I 2n — / . a o-0 I 2n • 



Proof. If {2n— 1, 2n} is a part of both odd-even pairings p and 7r, then let p' and 7r' denote the 
odd-even pairings obtained by deleting this part. Then (p,7r) g = q(p', 7r') 3 , so (p, LHS) g = 
(p, RHS) r Now suppose that {2n — 1, 2n} is a part of 7r but not p, and that instead 
{2n — l,w} and {6, 2n} are parts of p. Let p' denote the odd-even pairing obtained from p 
by deleting these two parts and replacing them with {b, to}, and let tt' be as above. Then 
(p', vr% = (p, ir) g , so (p, LHS> 9 = (p, RHS> 9 . □ 

Lemmas 13.71 and 13.91 imply that the left-hand-side and right-hand-side of transformation 

2 

Rule 2 are equivalent (=). 

Theorem 3.10. For any odd-even pairing t, there is an equivalent integer-linear combina- 

2 

tion of planar pairings r = a a a, with a a G Z. 

Proof. We start with a particular planar pairing, say tt ^= 2 1 4 1 6 I " " ' - ^ e nave the equivalence 

2 

71 = 71. We shall then make modifications on the left side, doing adjacent transpositions of 
even labels, making the same modification on the right, and then retransform the right side 
into a combination of planar pairings. Eventually we have converted the LHS into r and the 
RHS is an integer linear combination of planar pairings. 

Suppose we swap the labels b and 6 + 2. Assuming that the RHS was a linear combination 
of planar pairings, we need to check that after the swap the RHS can be transformed to again 
be a linear combination of planar pairings. Consider one such planar pairing of the RHS. 
If b + 1 were paired with one of b or b + 2, then after the swap the result is still planar. 
Otherwise the three chords containing b, 6 + 1, and 6 + 2 form three parallel crossings. These 
divide the remaining vertices into four contiguous regions, so that remaining chords only 
connect vertices in the same region. After the b,b + 2 swap, the three parallel crossings form 
an asterisk. The asterisk can be transformed using Rule 2 (Lemmas 13.71 and 13. 9p . None of 
the remaining chords cross the newly transformed chords, so once they are added back in, 
the results are planar. □ 



For an odd-even pairing r, let ^ CT a a o be any linear combination of planar pairings equiv- 
alent (=) to r. Now ^2 a (Vir^^ — a a )a lies in the null-space of A^2, but A^2 is nonsingular, 
so vjj?^ = a a for each a. In particular, the linear combination promised by lTheorem 3.101 is 
unique, and gives the r th column of p( DD ). This linear combination was obtained by repeated 
application of Rule 2, which completes the proof of lTheorem 1.41 which implies [Theorem 1.31 
as a corollary. 
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4 Comparing grove and double-dimer polynomials 

As we shall see in this section, the grove partition "L" polynomials and double-dimer pairing 
"X" polynomials are closely related to each other. In fact, the U X" polynomials are a 
specialization of the "L" polynomials. We start with a lemma about groves. 

Lemma 4.1. If a planar partition o contains only singleton and doubleton parts, and a' is 
the partition obtained from a by deleting all the singleton parts, then the "L" polynomials for 
a and a 1 are lexicographically equal. (We say "lexicographically equal" rather than "equal" 
because the formulas look the same even though the underlying "L variables" represent dif- 
ferent electrical quantities.) 

For example, Pr(13|2|4|5|68|7|9) = L li3 L 6>8 - L lt6 L 3 , 8 and Pr(13|68) = L^L 6fi -L lfi L^ 
but the Us mean different things in these two equations, since they are the current responses 
when different numbers of nodes are held at zero volts. Note that for general partitions we 
cannot simply drop the singleton parts. For example, Pr(123) = Li 2 iv2,3 + -^1,2-^1,3 + -^1,3-^2,3, 
while Pr(123|4) has these three terms together with a fourth term, namely Li^Z^. 

Proof. Recall the algorithm that computes columns of the projection matrix by finding 
for a given partition r an equivalent linear combination of planar partitions (the proof of 
ITheorem 2.6j) . Each partition in the result will have the same number of parts as r, and every 
singleton of r is a singleton of each partition in the combination. If we are only interested in 
the 0" th row of the projection matrix, where a has k parts of size 2 and rest singletons, then 
we need only consider columns r for which r has n — k parts. If r has any part of size 3 or 
more, then r has more singletons that a, and it does not contribute to cr's row. The only 
partitions r that contribute have k parts of size 2, and the same singletons that a does. For 
such partitions r, we may drop the singleton parts, find the equivalent linear combination 
of planar partitions, and then re-adjoin the singleton parts. Thus the computation of Pr(cr) 
precisely mirrors the computation of Pr(<r'). □ 

The following theorem is a reformulation of ITheorem 1.51 

Theorem 4.2. If a planar partition a contains only doubleton parts, and we make the 
following substitutions to the grove partition polynomial Pr (a): 

{0, if i and j have the same parity, 

(_l)(K-il-i)/2x ij3 ., otherwise, 

then the result is (— 1) CT times the double-dimer pairing polynomial Pr (a), when we interpret 
a as a pairing, and ( — 1)°" is the signature of the permutation a%, 03, ... , o^n-i- 

Proof. Consider the computation of Pr(cx) when the planar partition a contains only dou- 
bleton parts. When we express a partition r as a linear combination of planar partitions, 
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any singleton parts of r show up in each planar partition (with nonzero coefficient), so if r 
contains singleton parts, Va)r = 0. When we transform partitions according to Rule 1, the 
number of parts is conserved, so partitions that end up contributing to the a th row of the 
projection matrix will contain only doubleton parts. Rule 1 transforms a partition into a 
combination of six partitions, but the first four of them have the wrong part sizes to con- 
tribute to Va)r, so we may keep only keep the last two. In other words, we may use the 
abbreviated transformation rule 

13|24 — 12|34 — 14|23. 

Let us transform 14 1 25 1 36 using this abbreviated rule: 

14|25|36 -12|45|36 - 15|24|36 

-> -12|45|36 + 15|23|46 + 15|26|34 

-> -12|45|36 - 14|23|56 - 16|23|45 - 16|25|34 - 12|56|34. 
Compare this to the corresponding rule for double-dimer pairings: 

14|25|36 -> +12|45|36 + 14|23|56 - 16|23|45 + 16|25|34 - 12|56|34. 

The transformation rules are the same except for a sign factor equal to the signature of 
the pairing. Thus P£ D) = (-l) a 1t P$-. Recall now that the double-dimer projection 
matrix gave coefficients for monomials weighted by (— 1)* pairs that cross , and that this sign 
was by f|T2l equal to the signature of the pairing times the product over pairs in the 
pairing of (-lJdWI-iJ/z, □ 



5 Scaling limits and multichordal SLE 

The results in the preceding sections are all exact results that hold for any planar graph. 
In this section we study the asymptotic behavior of the partition probabilities when the 
graphs approximate a fine lattice restricted to a domain in the plane. We consider two 
different limits for groves; the first one gives multichordal loop-erased random walk, which 
in the limit converges to multichordal SLE2, while the second one in the limit converges 
to multichordal SLEg. We also obtain the limiting pairing probabilities of the paths in the 
double-dimer model, and the limiting connection probabilities for the contour lines of the 
discrete Gaussian free field with certain boundary conditions, and in particular show that 
they are equal. The latter of these two models converges to multichordal SLE 4 [SS06J [SS07J. 

Dubedat [Dub06j has also studied the pairing probabilities in multichordal SLE K . We only 
treat the 2, 4, 8, while Dubedat's calculations are relevant for the continuous range 

< k < 8. On the other hand, our calculations are more explicit. For example, Dubedat's 
solution is actually a multiparameter family of solutions that solve a certain PDE. (It is not 
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known a priori that the solutions in [Dub06] span the solution space, though for a given 
number of strands this can be verified a posteriori.) Each of these solutions is relevant to 
multichordal SLE K , but when there are more than two strands, it is not clear which solution 
to the PDE is the "canonical solution" that describes the behavior of a given discrete model 
such as spanning trees. Our calculations avoid that issue entirely by working directly with 
the discrete models. 



5.1 Multichordal loop-erased random walk and SLE2 

We consider multichordal SLE 2 on a domain U, with m chords connecting 2m points. To get 
SLE2 as a limit of groves, we consider groves on finer and finer planar grids approximating U, 
with n = 4m nodes, where each of the 2m SLE-endpoints is associated to a node, and each of 
the 2m intervals between the SLE-endpoints is wired together and associated to a node. (See 
Figure 8j ) The "interval" nodes are required to be in singleton parts of the corresponding 
grove partition, and the point nodes are required to be in doubleton parts. By lLemma 4.11 to 




Figure 8: Multichordal loop-erased random walk: Shown here is a random grove with 16 
nodes, conditioned to have the long extended (and unnumbered) nodes in singleton parts, 
and the short (numbered) nodes each in a part of size 2. In this case the pairing is \ \ \ | 1 1 g , 
and in the scaling limit for this geometry this pairing occurs with probability 9 4 6 8 5 7 7 7 7 7 6 — = 
0.0239. Within the grove the paths connecting the numbered nodes are shown in bold. Each 
of these connecting paths is, conditional upon the other paths, the loop-erasure of a random 
walk started at one endpoint and conditioned to hit the other endpoint before hitting either 
the boundary or the other paths. 
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compute the SLE 2 connection probabilities it is sufficient to compute the current responses 
amongst the point nodes (when the interval nodes are grounded), and substitute these values 
into the "L" polynomials associated with pairings of these 2m point nodes. 

By conformal invariance, we may without loss of generality fix U to be the upper half- 
plane, with point nodes at Xi, . . . , £ 2m € R, and we assume that no point node is at oo. 
We approximate the upper half-plane with the upper half cartesian lattice eL x eN, and 
round Xj to the nearest lattice point xf' G eTL x eN, and let each edge have unit resistance. 
The current response between distinct nodes xf and x^ is (assuming the nodes are non- 
adjacent) the voltage at the vertex one lattice spacing above x^ when x^ is at one volt and 
the remainder of the real axis is at zero volts, which is the probability that a random walk 
started one lattice spacing above Xj ends up at x^ when it first hits the real axis, which is 

= {I + o{l))- £ " (13) 

7T (Xi - XjY 



see e.g. |Spi76[ Chapter III §15], where the o(l) term goes to as £ goes to 0. 

When we consider the ratios of these probabilities, the e 2 /tt factors drop out and we can 
then take the limit e —>■ 0. We find then that Lij is inversely proportional to the square of 
the distance between points Xi and Xj. 

In the bichordal case, the normalized probabilities are 

Pr(12|34) =L 12 | 34 - L im 
Pr(14|23) =L 14 | 23 - L im . 

In the e — > limit, the unnormalized probability (conditional on there being two chords) is 
Pr(14|23) 



Pr(14|23) (xi— X4) 2 (X2—X3) 2 (X1—X3) 2 (X2—X4) 2 



Pr(14|23) + Pr(12|34) , 1 q < ,2 + ? V 2 ? ^ v 

vi/ 1 v 1 > (x 1 ~x 2 ) z (x 3 -x4y (xi-x4) z (x2-x 3 y (xi— X3Y {X2-X4,y 

If we let s denote the cross ratio s = (x 4 — x 3 )(x 2 — Xi)/[(x 4 — x 2 )(x 3 — xi)], this limiting 
probability can be written as 

2s 3 - s 4 



l-2s + 4s 3 -2s 4 ' 



which agrees with the known formula for bichordal SLE 2 , see [BBK05J § 8.2 and 8.3] [SW05] 
|Dub061 § 4.1 and 4.2]. 
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For trichordal SLE2, the relevant polynomials are 



Pr(12|34|56) 


— — -^15|26|34 


+ -^12|56|34 


+ -^14|26|35 + -^15|24|36 






— -^14|25|36 


+ -^13|25|46 


— -^12|35|46 — -^13|24|56 




Pr(16|23|45) 


= — -^16|24|35 


+ -^14|26|35 


+ -^15|24|36 — -^14|25|36 






+ -^16|23|45 


— -^13|26|45 


— -^15|23|46 + -^13|25|46 




Pr(14|23|56) 


=-^15|24|36 — 


-^14|25|36 — 


-^15|23|46 + -^13|25|46 + -^14|23|56 


— -^13|24|56 


Pr(25|16|34) 


=-^16|25|34 — 


-^15|26|34 — 


-^16|24|35 + -^14|26|35 + -^15|24|36 


— -^14|25|36 


Pr(36|12|45) 


=-^14|26|35 — 


-^12|46|35 — 


-^14|25|36 — -^13|26|45 + -^12|36|45 


+ -^13|25|46 



5.2 Grove Peano curves and multichordal SLEg 

We can do a similar calculation for multichordal SLEg containing m chords connecting 2m 
points. Lawler, Schramm, and Werner [LSW04] showed that on a fine grid in which part of 
the boundary is wired and part is free, the path separating a random spanning tree from its 
dual spanning tree converges to chordal SLE 8 connecting the two boundary points where the 
boundary conditions change from wired to free. If the boundary conditions alternate between 
wired and free and back m times, there will be m paths separating a random grove from its 




1 




Figure 9: The multichordal Peano curves (which converge to SLEg) associated to a grove 
with four nodes (the extended black segments). In this case the grove partition is 1|2|34, 
which gives the Peano curve pairing \ \ \ \ \ \ \ . In the scaling limit for this geometry these 
connections occur with probability (2 — y/2)/8 = 0.0732. 
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dual grove which connect the 2m points where the boundary conditions change from wired 
to free, and it is not difficult to see that each such path, conditional upon the other paths, 
converges to SLE 8 . Here the groves have n = m nodes, which correspond to "interval" nodes 
between the (2i — l) th and (2z) th endpoints of the SLE 8 strands. (See |Figure 9| ) We again 



consider U to be the upper half-plane and x%, . . . , x 2m the endpoints of the m SLE 8 strands. 
Let Ii be the interval node whose endpoints are x 2 i-\ and x 2 i. Now Ljj is the current into Ij 
for the harmonic function which is 1 on Ij, zero on the other ly and has Neumann boundary 
conditions elsewhere. 

For m — 2, the event that X\ . X4. cLFG connected and x 2 , £3 are connected is the event that 
the grove has one tree connecting nodes 1 (the interval I\ = [xx,x 2 ]) and 2 (the interval 
J 2 = [x 3 ,x 4 ]), which has probability Pr(12)/(Pr(12) + Pr(l|2)) = Li )2 /(L 1)2 + 1), where L lj2 
is the modulus of the rectangle which is the conformal image of U and whose vertices are 
the images of Xi, x 2 , 23, £4. By the Schwarz-Christoffel formula, the Riemann map from U 
to the rectangle is 

/(*) = 

Using f(z) we may express 



dw 



i(w — XiY^iw — x 2 ) l / 2 {w — x 3 y/ 2 (w — X4) 1 / 2 
f(x 2 ) - /(Xi) 



Ll ' 2= -i(f(xs)-f(x 2 )Y 
For m = 3, an interesting special case is where the domain is the regular hexagon and the 
three nodes are three nonadjacent sides of the hexagon. Recall now the dual electric network 



with its dual grove, which was illustrated in Figure 5 Here the dual electric network also 



converges to the regular hexagon, but with extended nodes on the other three nonadjacent 
sides of the hexagon. Thus the probability that the primal grove has type 1|2|3 equals the 
probability that the dual grove has type 1|2|3, but this latter event is the event that the primal 
grove has type 123, so Pr[l|2|3] = Pr[123]. From this we see 1 = Pr(l|2|3) = Pr(123) = 
£1,2-^2,3 + £1,2-^1,3 + £1,3-^2,3, and hence by symmetry Li j2 = Li j3 = L 2j3 = l/y/3. We can now 
compute the other normalized probabilities Pr(12|3) = Pr(13|2) = Pr(23|l) = L 2 ,3 = 1/VS- 
Thus the unnormalized partition probabilities are Pr[l|2|3] = Pr[123] = l/(2xl+3xl/v / 3) = 
2-^3 = 0.268 and Pr[12|3] = Pr[13|2] = Pr[l|23] = 2/^3-1 = 0.154. (The corresponding 
calculation for the connection probabilities for percolation (SLEg) on a regular hexagon was 
carried out by Dubedat [Dub06l § 4.4], but while the SLEg calculations here are almost trivial 
in retrospect, the corresponding SLE 6 calculations are much more intricate.) The scaling 



limit of the L's can similarly be worked out using duality for the example in Figure 9 
L iji+ i = 1/2 and L i i+2 = l/v2 — 1/2. The response matrix for the regular 2n-gon with sides 
that are alternately free and wired has a nice formula that is given in [BSW07J. 

When m > 2, the quantities Lij for more general domains can be represented geometri- 



cally as the moduli of images of maps of U onto various vertical slit rectangles, see Figure 10 
(In the remainder of this subsection i = We explain here how to compute the Life's 

when the domain U is the upper half-plane. Consider the nodes formed from the intervals 



Boundary partitions in trees and dimers 



Kenyon <fc Wilson 



35 



[xi, x 2 ], [x 3 , x 4 ], • • • , [x 2m -i, a?2m], where x t < x 2 < ■ ■ 
f defined in the upper half-plane by 



< x 2m . Consider the analytic function 



UT=i 2 ( w - a i> dw 



where a±, 



are parameters. By the Schwarz-Christoffel formula, if the a's are chosen 
judiciously, then / will map the upper half-plane to a rectangle with vertical slits, such that 
one of the intervals [x 2 j-i, %2j\ (j 1) is mapped to the top of the rectangle, the remaining 
intervals are mapped to the bottom of the rectangle, and the spaces between the intervals 
are mapped to the two sides of the rectangle and the vertical slits on the bottom of the 
rectangle. (The example in Figure 10 is for j = 2. The case j = 1 is similar, except that 
"top" and "bottom" are reversed.) The quantity Lj^ is then given by length of the image 
of [x2k-i,%2k\ divided by the height of the rectangle: 



f(x 2 k) ~ f{x2k-l 



(Here the choice of the a's, and hence /, depends upon j.) 
Define the analytic function g p in the upper half-plane by 

w p dw 



■ 



2m 



[W 



Each of the above analytic functions / = fj is a linear combination of go, gi, 
us define two m x (m — 1) matrices V and /' by 



, Qm-2- 



Let 



lm(g p (x 2 ) - g p {x 2 k-i)) and I' k = g p (x 2j ) - g p (x 2j 



2j-lj 




f{xn) /N fbi 



Figure 10: The function f(z) = £f ^i^L j- maps the upper half-plane to the slit 



rectangle, for suitable choices of the a/s. In the above figure L 



X3,X4\,[X7,X8\ 



'/(*s)-/(*a)" 
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Since for any j ^ 1, we can express fj as a linear combination of the g p 's, this same linear 
combination of the columns of V will be nonzero (the height of the slit rectangle) at row j 
and zero elsewhere, and the same linear combination of the columns of I' gives the lengths 
of the images of the intervals in the slit rectangle. Thus V has rank m — 1, and since the 
first row is all 0, we may adjoin an m th column which is all 1 to get a nonsingular matrix V, 
and adjoin an m th column to I' which is all to get a matrix /. Then for any vector v, we 
have LVv = Iv, which allows us to compute L by 

L = IV' 1 . 

Thus we may compute the limiting connection probabilities completely mechanically. 

5.3 Double-dimer crossings — mult ichor dal SLE4? 

There is a lot of evidence to support the hypothesis that the scaling limit of chains in the 
double-dimer model is given by SLE4, but this has not yet been rigorously proved. The 
formulas given here, which are for double-dimer path scaling limits, agree with the formulas 
for SLE 4 (see also ISection 5.41) . and so provide further evidence supporting this hypothesis. 

As in ISection 5.11 we approximate the upper half-plane with the upper half cartesian 
lattice eL x eN, and consider a set of black nodes near the points x\, x%, . . . , X2 m -i and white 
nodes near the points X2, x±, . . . , X2 m , letting node i be the nearest lattice point x^ G sZxeN 
to Xi which has the appropriate color. From |Ken00[ § 5.2] we have 

X;; = (l+ (1))-— ^— . (14) 

For the double-dimer model we can easily recover the cross-ratio formula that is known 
for bichordal SLE 4 |BBK05l § 8.2 and 8.3] [SW05] |Dub06l § 4.1] as follows: 

Pr(12|34) =X lj2 X 3 ,4 
Pr(14|23) =X M X 3 , 2 

we get 

1 1 

|xi-x 4 | IX2-X3I 

1 1_ 1 __L 1 

|xi-X2| |X4-X3| |xi-X4| |X2-X3| 

(x 2 - xi)(x 4 - x 3 ) 

(x 3 - X 1 )(x 4 ~ X 2 ) 

which is the formula for the cross-ratio. 



When x\ < x 2 < 23 < £4, in the scaling limit e 

Pr(14|23) = ^J^B 

Pr(12|34) + Pr(14|23) 
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For the trichordal case the relevant polynomials are 



Pr(12 


34 


56) 


=^1,2X3,4X5,6 


+ -^1,4-^3,6-^5,2 


Pr(16 


23 


45) 


=^1,6^3,2^5,4 


+ -^1,4-^3,6-^5,2 


Pr(14 


23 


56) 


=Xi,4^3,2^5,6 


— Xi^XsfiXsp 


P7(25 


16 


34) 


=-^1,6-^3,4-^5,2 


— Xi^Xs^Xsp 


Pr"(36 


12 


45) 


=X 1j2 X3 6 X 5j 4 


— X 1)4 X3 i6 X 5i2 



5.4 Gaussian free-field contour lines and multichordal SLE4 

Schramm and Sheffield studied the contour lines of the discrete Gaussian free field on a 
planar domain with boundary conditions +A on one segment of the boundary and —A on 
the rest of the boundary [SS06] . and they showed that for a certain value of A the contour 
line of height connecting the two special boundary points is in the scaling limit given 
by chordal SLE4. In general there can be multiple boundary points between which the 
boundary heights alternate between +A and —A, and the contour lines will form a random 
pairing of these boundary points. When there are multiple boundary points, at which the 
height alternates between +A and —A, each of these boundary points induces a drift term on 
the Brownian motion driving the SLE4 process |SS06j |SS07j . Given the conjecture that the 
scaling limit of the double-dimer model yields SLE 4 , and the fact that we know the scaling 
limit of the connection probabilities in the double-dimer model with multiple strands, it is 
natural to compare these probabilities to the connection probabilities of the contour lines 
in the discrete Gaussian free field. We already have a natural guess for these contour line 
connection probabilities (namely that they equal the connection probabilities in the double- 
dimer model), and it is not difficult to verify that this guess is in fact correct using the SLE 
formulation. In this subsection we prove 

Theorem 5.1. Consider the scaling limit of the discrete Gaussian free field on a planar 
lattice in the upper half-plane, with boundary heights that alternate between +A and —A at 
the points x\ < ■ ■ ■ < Xi n , for the value of A determined in JFS06J. Then the probability 
distribution of the manner in which the height-0 contour lines connect up with one another 
coincides with the connection probability distribution for the scaling limit of the double-dimer 
model in the upper half-plane with boundary nodes that alternate between black and white at 
locations xi, . . . , x 2n . 

Recall from lSection 3.ll and lSection 3.21 that the probability of any double-dimer pairing 
is a linear combination of Ds/D 's, where S ranges over balanced subsets of nodes. (An 
expression for D$ is recalled below.) We therefore start with a lemma giving the scaling 
limit of Ds/D . 
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Lemma 5.2. In the scaling limit of the double- dimer model in the upper half-plane, if the 
nodes are located at x± < ■ ■ ■ < X2 n o,nd the node colors alternate between black and white, 
and if S is a balanced subset of the nodes, then 



% = n i*i 



Proof. Recall that for a set S of nodes with equal numbers of odd and even nodes, we defined 
D s = det[(l i>i6S +l^ 5 ) x {-rf-^ 2 X h ^X l 



which, in the scaling limit, using Equation 14 when x\ < %2 < • • • < %2n, is given by 

j=l,3,...2n-l 



det 
±det 



[Ujes + Ujfts) x 



(_l)(|i-i|-l)/2 



i odd, i £ S 



det 



j even, j 6 S 



.7=2,4,.. .2n 
i odd, i £ S 

j even, j ^ S 



which is a product of two Cauchy determinants 



n ^ ~ 



i + j even 
i, j £ S or i,j ^ S 



n ^ ~ 

i + j odd 
i, j £ 5 or i, j £ S 

where the i,j G S factors are from the first determinant, and the i, j ^ S factors are from 
the second determinant. Since each multiplicand is positive and by lLemma 3.21 -Dq is a ratio 
of positive quantities, the global sign is +. 
We may evaluate the ratio 

TT~ = ~Y\_(- X i — x i) *' j 







where 



a 



if i,j G S or i,j <£ S, 



1,3 



if i G S and j ^ S or vice versa. 
We remark that products of the form nj<j( x j — x i)* i,j a l so arise in [Dub06j. 



□ 



Boundary partitions in trees and dimers 



Kenyon <fc Wilson 



39 



Proof of \ Theorem 5. 11 The strategy is to show that each D$/D is a martingale for the 
multichordal SLE4 diffusion associated with the contour lines of the discrete Gaussian free 
field | SS06 ] [SS07j . Then each of the double-dimer crossing probabilities will be a martingale 
for the diffusion, and since each of them has the right boundary conditions, this will show that 
the contour line crossing probabilities are equal to the double-dimer crossing probabilities. 
By symmetry considerations, it suffices to consider the SLE 4 process started at X\ where the 
other Xi's are force points. To check that M is a local martingale we need to verify that 

2 dM | dMy, (-1)'2 | 4 d 2 M 

1 ^ Xi — x\ dxi dxi ^—f Xi — Xi 2 dx 2 

The terms in the first sum come from the Loewner evolution of the points (i ^ 1) when 
the SLE 4 at X\ is developed, the terms in the second sum arise from the drifts that the 
points at Xi (i 7^ 1) induce on the Brownian motion driving the SLE4 at x±, and the last 
term comes from the quadratic variation of the Brownian motion driving the SLE4 process 
at x\. Upon dividing by 2M the equation becomes 

{) ^d t M _y> 1 d\ogM dhgMy^ d 2 hgM /dlogTWy 

2M ' Xi — X\ dxi dx\ ' Xi — x\ dx 2 V dx,\ J 

i>l i>l 1 

If M = Ds/D , then logM = Yl,i<j a i,j ^°s( x j ~ x i)i so the first term is 

(xk-xj^Xk-xx) i ^ k {x k -x j )(x j -x 1 ) i ^ k {x j -x l )(x k -x 1 ) ^{Xk-X^ 
the second term is 



the third term is 



2 ; 




«1,7 



and the fourth term is 



Adding these up we get 




aij-aij(-l) J - at! j + a ltj x ^ -a jjk - «ij(-l) fc - a lik (-l) j + 2a 1J a lfk 



Kj Kj<fc 



Xj-xi) 2 (xj - x 1 )(x k -an) 
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If aij 7^ 0, then ciij = (— l)- 7 , so the numerators in the first sum are all 0. For the second 
sum, there are a few cases to consider. If aj t t ^ then = (— and either a±j = 
or else Oi\^ = but not both, so the numerator simplifies to 0. If aij t k = then either 
a i,j — a i,k = (in which case the numerator simplifies to 0), or else a±j = (— l)- 7 and 
= (— and the numerator again simplifies to 0. Thus D s /D is a local martingale. 
Since D s /D is the probability that no node in S is paired to a node in S c (in the double- 
dimer model), it is bounded, and therefore it is actually a martingale. □ 
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A The response matrix 

To a pair (Q, N) consisting of a graph and a subset N of its vertices is associated an |N| x |N| 
matrix A, the response matrix, or Dirichlet-to- Neumann matrix. It is a real symmetric 
matrix defined from the Laplacian matrix as follows. The variables Lij used in this article 
are negatives of the entries of this matrix: L it j = — A it j. 



A.l Definition 

Let A be the Laplacian matrix of Q, defined on functions on the vertices of Q by A(f)(v) = 
^w^v c v,w(f(v) — f(w)), where the sum is over neighbors w of v and c v>w is the edge weight. 
Let N be the set of nodes of Q and I = V \ N the inner vertices. Then R v = R N © R 1 . In 
this ordering, we can write 

F 

H K. 



A = 



where K is an |I| x |I| matrix. 

Let V G R N be a choice of potentials at the nodes. If we hold the nodes at these 
potentials, the resulting current flow in the network will have divergence zero except at the 
nodes, that is, 

F G\ fv N \ /c N N 
H K) \ Vl ) [0 

where cn G R n is the vector of currents exiting the network at the nodes. 

In particular we must have v\ = —K~ 1 Hv-^ and so (F ' —GK~ l H)v-^ = cn- (The matrix K 
is invertible since it is the Dirichlet Laplacian on I with nonempty boundary N.) The matrix 
A — F — GK- l H satisfies 

Aw N = c N 
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and is called the response matrix, or Dirichlet-to-Neumann matrix of the pair (Q, N), and is 
also known as the Schur complement A/K. It is a symmetric positive semidefmite matrix 
from M N to M N with kernel consisting of the "all ones" vector [CdV98j . 

The following lemma is a special case of a more general determinant formula, but we 
need this special case in ISection 21 so we prove it here. 

Lemma A.l. Pr(tree)/ Pr(uncrossing) = detA, where A is obtained from A by removing 
any row and column. 

Proof. Since K is the Dirichlet Laplacian for the graph C7/N in which all vertices in N 
have been identified, by the matrix-tree theorem (often attributed to Kirchhoff, see [vLW92, 
Chapter 34]), det-fT is the weighted sum of spanning trees of G/~N. These lift to completely 
disconnected groves of Q, so det K = Z(uncrossing), where Z(t) denotes the weighted sum 
of groves of partition type r, in this case 1|2| • • • |n. 

Let A e = A+£/n; recall that this is the Dirichlet Laplacian for the graph Q' obtained from 
Q by adjoining an external vertex connected to every vertex in N by an edge of weight e, 
i.e., A e is obtained from A(Q') by striking out the row and column corresponding to the 
adjoined external vertex, so that lim e ^ det(A + e/ N )/e = nZ(tree). 

Let A e = A £ /K = F + eI N - GK^H. Then 

a _ (K GK-*\ // N 0\ 
* £ -{0 h J\H Kj> 

so det A £ = det A £ det.K" (a standard fact about Schur complements). 

But if A denotes A with a row and column removed, then det A = - lim e ^o(det A e )/e = 
^ lim e ^ (det A e )/e/ det K = Z (tree) jZ (uncrossing). □ 



A. 2 Relation between the L^j's and i^j's 

The A matrix is related to the matrix of pairwise resistances between the nodes as follows. 
We assume here for convenience that the electrical network is connected. 

Proposition A. 2. We have 

R id = (5 i -8 j ) T A- 1 (5 i -5 j ). 

Note that the right hand side is well defined since Si — Sj is in the image of A, and the 
inverse image of this vector is well-defined up to addition of an element of the kernel of A, 
which is perpendicular to Si — Sj. 

Proof. Find v € MX such that Av = —Si + Sj. This represents the (unique up to an additive 
constant) choice of potentials for which there is one unit of current flowing into the circuit 
at i and out at j. By definition, the resistance between i and j is the difference in potentials 
at i and j for this current flow. □ 
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The equation of Proposition A. 2 can be given in a computationally simpler form as 
follows. When finding v € M. v such that Av = —5i + Sj, we may choose v so that v n = 0, 
and let v denote v with v n dropped. Let A be the matrix obtained from A by deleting the 
last row and column. Then Av = —5i + 5j, where 5 n is interpreted as 0. By the matrix-tree 
theorem, A is invertible since det A is the weighted sum of spanning trees with edge weights 
given by — A i; j. Thus R itTl = A^j, and if i, j ^ n then 



Ri,j = (^i - 5j) T A 1 (8 i - 5j 



(15) 



Reciprocally to proposition A.2| we can write A in terms of R as follows: 
Proposition A. 3. We have A~l = Ri jU , and 



A; 



l :J 



r(-Ri,n Rj,n -Ri,i)" 



Proof. This is just the inverse of the transformation in ([15 



□ 



B Complete list of grove probabilities on 4 nodes 

Here we give the formulas for the planar partition probabilities when there are 4 nodes. 
There are 14 planar partitions on 4 nodes, but for each class of planar partitions that may 
be obtained from one another by cyclically rotating and/or reversing the indices, we need 
only give the formula for one representative planar partition in the class. 



P?(1234) =1, 

P7(2|143) =1^2,3 + 1^1,2-^1,3, 
P7(14|23) =1^,3 + |i?2,4 - ~#1,4 - |^2,3, 

Pr(l|2|34) =(+-Rl,2-Rl,3 + -^1,2-^2,3 + -Rl,2-Rl,4 + Rl, 2-^2,4 ~ 2^1,2-^3,4 ~ -^1,2 

+ -^1,3-^2,4 + -^1,4-^2,3 — -Rl,3-Rl,4 _ -^2,3-^2,4) /4, 
Pr(l|3|24) =(+-Rl,3-Rl,2 + -^1,3-^2,3 + -Rl,3-Rl,4 + -Rl, 3-^3,4 ~ 2^1,3^2,4 ~ R\^ 

+ -Rl, 2-^3,4 + -Rl,4-R2,3 — -Rl,2-Rl,4 ~ R2^R?,,a) /4, 
Pr(l|2|3|4) =(+-Rl,2-R2,3-R3,4 + -Rl,2-R2,4-R4,3 + -^1,3-^3,2-^2,4 + -^1,3-^3,4-^4,2 

+ -#1,4-^4,2-^2,3 + -Rl,4-R4,3-R3,2 + -R2,l-Rl,3-R3,4 + -R2,l-Rl,4-R4,3 

+ -R2,3-R3,l-Rl,4 + -R2,4-R4,l-Rl,3 + -#3,1^1,2 -#2,4 + -R3,2-R2,l-Rl,4 

— -Rl,2-R2,3-R3,l — -Rl,2-R2,4-R4,l — -Rl,3-R3,4-R4,l — -R2,3-R3,4-R4,2 

— Rl^R-s^ — -#l t 2-#3,4 — -#l,3-#2,4 — -#l,3-#2,4 — -#l,4-#2,3 — -#l,4-#2,3) /4. 
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Equivalent ly, 



Pr(1234) 




-^1,2-^2,3-^3,4 + -^1,3-^2,3-^2,4 + -^1,2-^2,4-^3,4 + -^1,4-^2,3-^2,4+ 
-^1,2-^1,3-^2,4 + -^1,2-^1,4-^2,3 + -^1,3-^2,4-^3,4 + -^1,4-^2,3-^3,4 + 
-^1,3-^1,4-^2,3 + -^1,2-^1,3-^3,4 + -^1,3-^1,4-^2,4 + -^1,2-^1,4-^3,4, 



Pr(2|143) 
Pr(14|23) 
Pr(l|2|34) 
Pr(l|3|24) 
Pr(l|2|3|4) 



-^2|143 + -^13|24 =-^1,3-^1,4 + -^1,4-^3,4 + -^1,3-^3,4 + -^1,3-^2,4, 

-^14|23 — -^13|24 =-^1,4-^2,3 — -^1,3-^2,4, 

-^1|2|34 =-^3,4, 

-^1|3|24 =-^2,4, 

-^1|2|3|4 =1- 
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